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Abstract We propose to interpret the stellar IMF as a property of the turbulence 
in the star-forming gas. Gravitationally unstable density enhancements 
in the turbulent flow collapse and form stars. Their mass distribution 
can be derived analytically from the power spectrum of the turbulent 
flow and the isothermal shock jump conditions in the magnetized gas. 
For a power spectrum index f3 = 1.74, consistent with Larson's veloc- 
ity dispersion-size relation as well as with new numerical and analytic 
results on supersonic turbulence, we obtain a power law mass distribu- 
tion of dense cores with a slope equal to 3/(4 — 0) = 1.33, consistent 
with the slope of Salpeter's stellar IMF. Below one solar mass, the mass 
distribution flattens and turns around at a fraction of a solar mass, as 
observed for the stellar IMF in a number of stellar clusters, because 
only the densest cores are gravitationally unstable. The mass distribu- 
tion at low masses is determined by the Log-Normal distribution of the 
gas density. The intermittent nature of this distribution is responsible 
for the generation of a significant number of collapsing cores of brown 
dwarf mass. 



1. Introduction 

The origin of the stellar initial mass function (IMF) is a fundamental 
problem in astrophysics because the stellar IMF determines photometric 
properties of galaxies and the dynamical and chemical evolutions of their 
interstellar medium. In this contribution we address the relation between 
statistical properties of turbulence and the origin of the stellar IMF as 
discussed in Padoan & Nordlund (2002, 2004). The main result of these 
works is that the power law slope, s, of the stellar IMF measured by 
Salpeter (1955) is the consequence of the turbulent nature of the star- 
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forming gas and is directly related to the turbulent power spectrum 
slope, (3, s = 3/(4 — /?). From this point of view it is not surprising 
that the origin of Salpeter's result has remained mysterious for half a 
century, as our understanding of turbulence has not improved much since 
the seminal work by Kolmogorov (1941). The situation has changed 
during the last decade, because ever increasing computer resources have 
recently allowed significant progress in both fields of turbulence and star 
formation. 

Since Salpeter's work, the stellar IMF has been measured successfully 
in many systems. We now know that the IMF in young clusters (e.g. 
Chabrier 2003) reaches a maximum at a fraction of a solar mass, and 
then turns around with a relative abundance of brown dwarfs (BDs) that 
may vary from cluster to cluster (e.g. Luhman et al. 2000). The work 
we present here explains also this feature of the IMF as a consequence 
of the supersonic turbulence in the star-forming gas. 

Using the properties of supersonic turbulence we have derived from 
recent numerical simulations, we predict the stellar IMF essentially with- 
out free parameters. This predicted IMF is shown to depend on the rms 
Mach number, mean density and temperature of the turbulent flow. It 
is also shown to agree well with Salpeter's IMF for large stellar masses 
and with the low mass IMF derived for young stellar clusters. Our view 
of the IMF as a natural property of supersonic turbulence in the mag- 
netized and isothermal star-forming gas provides an explanation for the 
origin of BDs as well. According to this picture, BDs may be formed in 
the same way as hydrogen-burning stars. 

2. Statistics of Supersonic MHD Turbulence 

The velocity power spectrum in the inertial range of turbulence (be- 
tween the scales of energy injection and viscous dissipation) is a power 
law, E v (k) oc k~P , where k is the wave-number. The spectral index is 
(3 « 5/3 for incompressible turbulence (Kolmogorov 1941), and j3 « 2 for 
pressureless turbulence (Burgers 1974; Gotoh & Kraichnanl993). In re- 
cent numerical simulations of isothermal, super-Alfvenic and highly su- 
personic magneto-hydrodynamic (MHD) turbulence we have obtained a 
power spectrum intermediate between the Burgers and the Kolmogorov 
power spectra (Boldyrev, Nordlund & Padoan 2002) and consistent with 
the prediction by Boldyrev (2002), E(k) oc A;" 1 - 74 . 

The probability density function (PDF) of gas density in isothermal 
turbulent flows is well approximated by a Log-Normal distribution with 
moments depending on the rms Mach number of the flow (Nordlund &; 
Padoan 1999; Ostriker, Gammie & Stone 1999). The density structure 
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is characterized by a complex system of interacting shocks resulting in a 
fractal network of dense cores, filaments, sheets and low density "voids" , 
with a large density contrast. Most of the mass concentrates in a small 
volume fraction (according to the Log-Normal PDF), a manifestation 
of the intermittent nature of the turbulence. An example of a projected 
turbulent density field is shown in Figure 1. 




Figure 1. Projected density field from a numerical simulation of isothermal, super- 
sonic hydrodynamic turbulence with an effective resolution of 1024 3 computational 
zones (Kritsuk, Padoan & Norman, in preparation). The contrast has been reduced 
in order to show details of the low density regions. 

3. From Kolmogorov to Salpeter: The Mass 
Distribution of Collapsing Cores 

A simple model of the expected mass distribution of dense cores gener- 
ated by supersonic turbulence has been proposed in Padoan & Nordlund 
(2002). The model is based on two assumptions: i) The power spectrum 
of the turbulence is a power law; ii) the typical size of a dense core 
scales as the thickness of the postshock gas. The first assumption is a 
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basic result for turbulent flows and holds also in the supersonic regime 
(Boldyrev et al. 2002). The second assumption is suggested by the fact 
that postshock condensations are assembled by the turbulent flow in a 
dynamical time. Condensations of virtually any size can therefore be 
formed, independent of their Jeans' mass. 

With these assumptions, together with the jump conditions for MHD 
shocks, the mass distribution of dense cores can be related to the power 
spectrum of turbulent velocity, E v (k) oc k' 13 : 



If the turbulence spectral index (3 is taken from the analytical predic- 
tion by Boldyrev (2002), which is consistent with the observed velocity 
dispersion-size Larson relation (Larson 1979, 1981) and with our numer- 
ical results (Boldyrev et al. 2002), then « 1.74 and the mass distribu- 
tion is N(m) dlnm oc m.~ L33 dlnm, almost identical to Salpeter's stellar 
IMF (Salpeter 1955). The exponent of the mass distribution is rather 
well constrained, because the value of (3 for supersonic turbulence can- 
not be smaller than the incompressible value, (3 = 1.67 (slightly larger 
with intermittency corrections), and the Burgers case, (3 = 2.0. As a 
result, the exponent of the mass distribution is predicted to be within 
the range of values of 1.3 and 1.5. 

While massive cores are usually larger than their critical Bonnor- 
Ebert mass, tobe, the probability that small cores are dense enough to 
collapse is determined by the statistical distribution of core density. In 
order to compute this collapse probability for small cores, we assume i) 
the distribution of core density can be approximated by the Log-Normal 
PDF of gas density and ii) the core density and mass are statistically 
independent. Because of the intermittent nature of the Log-Normal 
PDF, even very small (sub-stellar) cores have a finite chance to be dense 
enough to collapse. Based on the first assumption, we can compute the 
distribution of the critical mass, p(?71be) drriBE, from the Log-Normal 
PDF of gas density assuming constant temperature (Padoan, Nordlund 
& Jones 1997). The fraction of cores of mass m larger than their critical 
mass is given by the integral of p(t«be) from to m. Using the second 
assumption of statistical independence of core density and mass, the 
mass distribution of collapsing cores is 



This mass distribution is a function of the rms Mach number, mean 
density and temperature of the turbulent flow, as these parameters enter 
the PDF of gas density and thus p{m^). Figure 2 (left panel) shows 



N(m) dlnm cc m ^dlnm . 



(1) 
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Figure 2. Left panel: Analytical mass distributions computed for (n) = 10 4 cm~ 3 , 
T — 10 K and for three values of the sonic rms Mach number, Ms = 5, 10 and 20 
(solid lines). The dotted lines show the mass distribution for T = 10 K, Ms = 10 and 
(n) = 5 x 10 3 cm -3 (lower plot) and (n) = 2 x 10 4 cm -3 (upper plot). Right panel: 
IMF of the cluster IC 348 in Perseus obtained by Luhman et al. (2003) (solid line 
histogram) and theoretical IMF computed for (n) = 5 x 10 4 cm~ 3 , T — 10 K and 
M s = 7 (dashed line). 



five mass distributions computed from equation (2) with three different 
values of the sonic rms Mach number and two different values of density. 
Padoan & Nordlund (2004) have suggested that BDs may originate from 
the process of turbulent fragmentation like hydrogen-burning stars. This 
is illustrated by the analytical IMF in the left panel of Figure 2, which 
shows a relatively large abundance of BDs is predicted for sufficiently 
large values of the mean density or rms Mach number. 

The IMF of the cluster IC 348 in Perseus, obtained by Luhman et al. 
(2003), is plotted in the right panel of Figure 2 (solid line histogram). 
The IMF of this cluster has been chosen for the comparison with the 
theoretical model because it is probably the most reliable observational 
IMF including both BDs and hydrogen burning stars. In Figure 2 
we have also plotted the theoretical mass distribution computed for 
(n) = 5 x 10 4 cm" 3 , T = 10 K and Ms = 7. These parameters are appro- 
priate for the central 5x5 arcmin of the cluster (0.35 x 0.35 pc), where 
the stellar density corresponds to approximately 2 x 10 4 cm~ 3 . The fig- 
ure shows that the theoretical distribution of collapsing cores, computed 
with parameters inferred from the observational data, is roughly consis- 
tent with the observed stellar IMF in the cluster IC 348. Similar IMFs 
have been obtained for several other young clusters (Chabrier 2003). 

4. Conclusions 

We have proposed to explain the stellar IMF as a property of super- 
sonic turbulence. This scenario is very different from previous theories of 
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star formation (see Shu, Adams & Lizano 1987), where it is assumed that 
stars of small and intermediate mass are formed from sub-critical cores 
evolving quasi-statically, on the time-scale of ambipolar drift. These 
theories do not account for the ubiquity of turbulence in star-forming 
clouds and therefore ignore the effect of turbulence in the fragmentation 
process. 

A naive interpretation of our results may lead to the conclusion that 
the stellar IMF at large masses should be a universal power law, with 
a slope very close to Salpeter's value. The statistics of turbulence con- 
trolling the origin of the stellar IMF are indeed universal (they depend 
on flow properties such as the rms Mach number, not fluid properties, 
and are insensitive to initial conditions that are soon "forgotten" due 
to the chaotic nature of the turbulence). However, such statistics are 
derived from ensemble averages. According to our derivation, massive 
stars originate from shocks on relatively large scales. In any given system 
(for example a molecular cloud) the number of large scale compressions 
(thus the number of massive stars) is relatively small. Because of the 
small size of the statistical sample and because of the intermittent na- 
ture of the turbulence, large deviations from Salpeter's IMF are possible 
in individual systems. Salpeter's IMF is therefore a universal property 
of supersonic turbulence, but it is not necessarily reproduced precisely 
in every star-forming region. 
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